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ABSTRACT 


A unidirectional fiberglass epoxy composite specimen is 
modelled as a homogeneous transversely isotropic continuum plate 
medium. Acousto-ultrasonic non-contact input-output characteri- 
zation by tracing SV waves in the continuum is studied theoreti- 
cally with a transmitting and a receiving transducer located on 
the same face of the plate. The isotropic plane of the equi- 
valent continuum plate model lies in the midplane of the plate 
and is parallel to the top and bottom faces of the plate. 

The single reflection problem for an incident SV wave at a 
stress-free plane boundary in a semi- infinite transversely iso- 
tropic medium whose isotropic plane is parallel to the plane 
boundary is analyzed first. For all angles of incidence, the 
angle of reflection of the SV wave is equal to the angle of inci- 
dence of the SV wave. It is found that an obliquely incident SV 
wave results in a reflected SV wave and a reflected P wave for an 
angle of incidence of the incident SV wave less than the critical 
angle 52°. when the angle of incidence of an incident SV wave is 
equal to or greater than the critical angle, there exists only an 
SV wave in the medium as the reflected P wave degenerates into a 
surface wave travelling parallel to the plane boundary. The 
amplitude of the surface wave decays exponentially with the per- 
pendicular distance from the plane boundary. The amplitude 


1 



ratios of the reflected P and SV waves to the incident SV wave as 


a function of the angle of incidence are plotted from zero to the 
critical angle. The amplitude ratio of the reflected SV wave is 
found to be minus one when the angle of incidence is equal to or 
greater than the critical angle. The balance in energy flux nor- 
mal to the plane boundary is checked. 

Accordingly, the delay time for wave propagation between the 
transmitting and the receiving transducers is computed as if the 
SV waves were propagating in a half space. It is found that the 
directional dependence of the phase velocity of the SV wave pro- 
pagating in the transversely isotropic medium has a significant 
effect on the delay time, as opposed to the directional indepen- 
dence of the phase velocity of a shear wave propagating in an 
isotropic medium. 

The displacements associated with the SV wave in the plate 
and which may be detected by the non-contact receiving transduc- 
er are approximated by an asymptotic solution for an infinite 
transversely isotropic medium subjected to a harmonic point load. 
The polar diagrams for the directivity function of the stresses 
due to SV waves in the plate are shown at frequencies of 0.75, 
1.50 and 2.25 MHz. 
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This study enhances the quantitative understanding of acous- 
to-ultrasonic non-destructive evaluation (NDE) parameters such as 
the stress wave factor (SWF) and wave propagation in fiber rein- 
forced polymeric, ceramic or metallic composites, which can be 
modelled as transversely isotropic media. 


3 



INTRODUCTION 


Fiber reinforced composite materials are attractive materials 
for aerospace applications because of their high specific 
mechanical properties. It has been shown that many composites, 
such as fiberglass epoxy composites or fiber reinforced ceramics 
as shown in Fig. 1, may be modelled as a homogeneous transversely 
isotropic continuum [1]. In this work, acousto-ultrasonic (AU) 
non-contact input-output characterization of a homogeneous trans- 
versely isotropic elastic plate is investigated by tracing SV 
waves . 

First, the single reflection problem of an incident SV wave 
at a stress -free plane boundary in a semi- inf inite transversely 
isotropic medium whose isotropic plane is parallel to the plane 
boundary is considered. At such boundaries, the conditions for 
the existence of wave mode conversion, critical angle phenomena, 
reflected surface wave, the angles of reflection of the reflected 
waves, and the amplitude ratios of the reflected waves to the 
incident wave are derived. 

Second, the SV wave input -output relations are derived when 
multiple reflections occur at the top and the bottom faces of the 
plate. The delay time between input and output versus the dis- 
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tance separating the transmitting and receiving transducers is 
analyzed. The directivity functions of the stresses associated 
with the SV waves are computed. And, the output displacement at 
the non- contact receiving transducer is approximated by an 
asymptotic solution. 

This investigation should enhance the quantitative understan- 
ding of AU NDE parameters such as the stress wave factor. It 
also provides the potential for assisting in the development of 
better NDE schemes utilizing the SWF as well as other AU parame- 
ters for the characterization of fiber reinforced polymeric, 
metal and ceramic composites. 
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SINGLE REFLECTION PROBLEM AT STRESS -FREE PLANE BOUNDARY IN 


SEMI-INFINITE TRANSVERSELY ISOTROPIC MEDIUM WHOSE ISOTROPIC PLANE 
IS PARALLEL TO PLANE BOUNDARY FOR INCIDENT SV WAVE 
1 . Reflected P and SV Waves 

For a homogeneous, linearly elastic transversely isotropic 
continuum, the number of independent elastic constants is five 
[1] . Define a coordinate system (x, y, z) for a semi- inf inite 
transversely isotropic medium whose isotropic plane is parallel 
to the plane boundary where the reflection occurs as follows: the 
plane boundary contains the x and y axes, and the z axis is the 
zonal axis of the medium, which is in the direction parallel to 
the fiber direction shown in Fig. 1. Also, see Fig. 2. The gen- 
eralized Hooke's law is written, relative to the (x, y, z) coor- 
dinate system, as [1] 

r xx " c ll u - x + c 12 v -y + c 13 w - z 

Tyy - C u n,x + C^v.y + C^w.z 

T zz - c 13 u - x + c 13 v -y + c 33 w > z 

r xz “ c 44( u - z + w,x) (1) 

r yz — C^^(v,z + w,y) 

r xy ” c 66< u -y + v » x ) 


where r rs (r, s - x, y and z) are the normal (r - s) and shear 
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(r * s) stresses with respect to the coordinate system (x, y, z); 
u, v and w are the displacement components of a point in the 
medium along the x, y and z axes, respectively; denotes par- 

tial differentiation with respect to the variable which follows; 
and C^i, C^ 2 > ^13 > C 33 anc * C 44 are five independent elastic 
constants where Cgg - 1/2(C]^ - C^ 2 ) . 

Let a progressive wave be represented as [2] 

(u, v, w) - A (P x , P y , P z ) exp{iu(S x x + S y y + S z z - t) } (2) 

where S x , Sy and S z are the components of the slowness vector, 
which is in the same direction as the normal to the wavefront and 
whose magnitude is equal to the reciprocal of the magnitude of 
the phase velocity [1], along the x, y and z axes, respectively; 
P x , Py and P z are the components of a unit vector of particle 
displacement along the x, y and z axes, respectively; A is the 
amplitude of particle displacement; t denotes time and u> denotes 
radian frequency. It follows from Eqs . (1) and (2) that the 
stresses can be represented as 

r xx - iwA[C u S x P x + C 12 s y p y + C 13 S z P z ]exp{iw(S x x + S y y + S z z - t) ) 

r yy - iwA[C 12 S x P x + C n S y Py + C 1 3 S z P z ] exp { 1 w ( S x x + s y y + S z z - t) } 

r zz = iwA[C 13 S x P x + C 13 S y Py + C 33 S z P z ]exp{iw(S x x + S y y + S z z - t) } 
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r xz - iwA[C 44 S z P x + C 44 S x P z ]exp{ iw(S x x + S y y + S z z - t) } 


(3) 


T yz - iwA[C44S z Py + C44S y P z ]exp{ iw(S x x + S y y + S z z - t) } 
r xy ” iu>A[CggSyP x + C^gS x Py ] exp { iw( S x x + S y y + S z z - t) ) 


The stress boundary conditions on the stress -free plane 
boundary require that [2] 


. (I) 

+ 

r (R> - 0 


xz 


1 xz 


(I) 

yz 

+ 

, yz ( R ) - 0 

(4) 

(I) 

zz 

+ 

rz m . 0 



where t rz ^ (r - x, y and z) represent stresses on the plane 
boundary associated with the incident SV wave, and r rz ^^ (r = x, 
y and z) represent stresses on the plane boundary associated with 
the reflected waves . 


In order to satisfy Eq. (4), it is required [2] that the 
frequency, u>, of the reflected wave be equal to that of the 
incident wave and that [2] 


S X (I) - s x W 


S (!) 


s < R ) 

y 


(5) 
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As a result of Eq. (5), the slowness vectors of the incident and 
reflected waves lie in a plane called the plane of incidence. 

This analysis can be simplified by assuming that the plane of 
incidence is the x — 0 plane; that is, the slowness vectors of 
the incident and reflected waves are in the x - 0 plane, as shown 
in Fig. 2. Then it follows from Eqs. (5) that 


s x (I) - S X (*> - o 


( 6 ) 


It has been shown [3] that, except along the principal mate- 
rial axes, P waves or SV waves travelling in a plane containing 
the zonal axis, z axis, of a transversely isotropic medium are 
quasi -longitudinal and quasi- transverse , respectively; that is, 
in general, the components of the unit vector of particle dis- 
placement for either the P or SV waves along the y and z axes, 

Py and P z , do not vanish; whereas the components along the x 
axis, P x , do vanish. Therefore, it follows from Eqs. (3) and (6) 
that stresses associated with P and SV waves are 


' xz 


xy 


r xx * 


0; r yy * 0; 


r zz * 


0 ; r yz * 0 


(7) 


It has also been shown that an SH wave travelling in a plane 
containing the zonal axis.z axis, of a transversely isotropic 
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medium possesses a transverse displacement only, that is, for the 
coordinates in Fig. 2, (P x , Py, P z ) - (1, 0, 0) [3]. Therefore, 
it follows from Eqs . (3) and (6) that the stresses associated 
with SH waves are 


r yz = r xx “ r yy “ r zz “ ® 

( 8 ) 

r xy " °! r xz * 0 

Assume that an SV wave is incident on the plane boundary, the 
x-y plane in Fig. 2. It follows from Eqs. (4) and (7) that 


r yz 


(I) 


* 0 ; 
r xz 


'zz (I) - 0 

(I) - 0 


(9) 


As a result of Eq. (9), it is known from Eq. (4) that t xz ( r ) is 
equal to zero. This means that no SH wave will be reflected back 
into the medium because a reflected wave of the SH type would 
result in a nonzero value of the stress r xz . So, Eqs. (4) 
reduce to 


, yz a> ♦ , y2 w - o 
+ ,„<» . o 


( 10 ) 
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Since either a reflected P wave or a reflected SV wave results in 
nonzero values of the stresses Ty Z ^ and r Z2 W, it is therefore 
concluded from Eqs . (7) and (10) that both a P wave and an SV 
wave may be reflected back into the medium. 

2 . Slowness Surface for P and SV Waves 

The equations of motion relative to the coordinate system (x, 
y, z) are [1] 

r xx, x + r xy,y + T xz,z “ £ u > tt: 

r xy , x + r yy ,y + T yz,z “ P v . tc ( n ) 

r xz, x + r yz,y + T zz,z ” 

where the body forces are identically zero for the homogeneous 
solution. 

It follows from Eqs. (1), (2) and (11) that the following 
equations of motion are obtained: 

[Cll s x 2 + c 66 s y 2 + c 44 s z 2 ' ^l p x + ( c 12 + c 66> s x s y p y 
+ (C]_3 + C 44 )S X S Z P Z ” 0 
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( 12 ) 


(Cl2 + Cgg)S X SyPjj + [CggS X ^ + C-^^Sy^ + - P ] Py 

+ (C]_3 + C 44 )SyS z P z “ ® 

(C 13 + C 44 )S X S Z P X + (C 13 + C 44 )S y S z Py + [C 44 (S X 2 + Sy 2 ) 

+ c 33 s z 2 -p]P z - 0 

The condition for the existence of the plane wave solution is 
expressed by setting the determinant of the matrix of the coef- 
ficients of P x , Py and P z in Eq. (12) equal to zero [1] : 


t c 11 s x 2 + c 66 s y 2 + c 44 s z 2 

(C 12 + c 66 )s x s y 

<c 13 + C 44 )S x S z 

(c 12 + c 66 )s x s y 

t c 66 s x 2 + c 11 s y 2 

+ ^44®z 2 " ^13 * ^44^®y®z 

(c 13 + c 44 )s x s z 

< c 13 + C 44 )S y S z 

tC 44 CS x 2 ♦ S y 2 ) ♦ c 33 s z 2 -p] 



(13) 

By expanding Eq. (13), 

three sheets 

of slowness surface are 


obtained. The slowness surface for an SV wave is given in [3] as 


C 11 + C 44 


(- 


)(S X 2 + S y 2 ) + ( 


c 44 + c 33 


•) s 5 


1 / 2 ' 


(C n - C 44 )(S X 2 + Sy 2 ) + (C 33 - C 44 )S Z 


(14) 


+ 4(S X 2 + S y 2 )S z 2 


( c ll ' c 44^ ( c 33 * c 44> 


1 1/2 


- (C ]_ 3 + c 44 )- 
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where S x , S y and S z are the components of the slowness vector for 
an SV wave along the x, y and z axes, respectively; p is density; 
and C]^, C^ 3 , C 33 and C 44 are elastic constants. Similarly, the 
slowness surface for a P wave is [3] 


C U + Caa 

( — — )(S X Z + S z ) + ( 

2 x y 2 


2 + Sy 2 ) + ( £ 33 _ti4 f L) s i 2 

2 a. c 2\ , r r - - _ r, . \ c 2 


+ l/2|^(C n - C 44 )(S X ^ + S/) +(C 33 - C 44 )S Z 
+4(S X 2 + S y 2 )S z 2 |(C n - C 44 ) (C 33 + C 44 ) 

' ( c 13 + c 44> 2 J| " P 


(15) 


where S x , S y and S z are the components of the slowness vector for 
a P wave along the x, y and z axes, respectively. 


One quadrant of the intersection of the slowness surface of 
an SV wave in the unidirectional fiberglass epoxy composite shown 
in Fig. 1 with the plane x - 0 and one quadrant of the intersec- 
tion of the slowness surface of a P wave in the unidirectional 
fiberglass epoxy composite with the plane x — 0 are shown from 
Eqs . (14) and (15) in Fig. 3; where the numerical values of the 
elastic constants and density given in [ 1 ] for the unidirectional 
fiberglass epoxy composite are used and are as follow: C-q — 
10.581 x 10 9 N/m 2 , C 13 - 4.679 x 10 9 N/m 2 , C 33 - 40.741 x 10 9 N/m 2 , 
C 44 - 4.422 x 10 9 N/ra 2 , and p - 1850 kg/m 3 . 
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3 . Angle of Reflection 


It follows from Eqs . (5) and (6) that the y-component of the 
slowness vector of an incident SV wave is equal to the y-compo- 
nent of the reflected SV wave as well as the y-component of the 
reflected P wave. Accordingly, the relation between the y-com- 
ponent of the slowness vector of an incident SV wave and the 
y-component of the slowness vectors of the reflected SV and P 
waves is given as 

Syd) - S y (SV) - Sy CP) - b (16) 

where Sy^ 1 ) represents the y-component of the slowness vector of 
an incident SV wave; S y ( sv ) represents the y-component of the 
slowness vector of the reflected SV wave; Sy( p ) represents the 
y-component of the slowness vector of the reflected P wave; and b 
is a common constant, as shown in Fig. 3. 

It follows from Eqs. (6), (14) and (16) that the relation 
between the z-component of the slowness vector of an incident SV 
wave and that of the reflected SV wave is 

S^ 1 ) - - S z < sv > (17) 

The minus sign is due to the fact that the slowness vector of an 
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incident SV wave points out of the medium, whereas the slowness 
vector of the reflected SV wave points into the medium, as shown 
in Fig. 2. Consequently, the value of the z- component of the 
slowness vector of an incident SV wave, S z ^^, is negative, 
whereas that of the reflected SV wave, S Z ^V), i s positive. 

The angle of reflection is defined as the angle between the 
slowness vector of a reflected wave, either type P or SV, and the 
normal to the plane boundary where the reflection occurs. Simi- 
larly, the angle of incidence is defined as the angle between the 
slowness vector of an incident SV wave and the normal to the 
plane boundary, as shown in Fig. 2. Therefore, the angle of 
reflection of a reflected SV wave is defined as 



and the angle of incidence of an incident SV wave is defined 
as 

6 Z - tan- 1 (S y < I >/-S z < I >) (19) 

It follows from Eqs. (16) through (19) that, for all angles, the 
angle of incidence of an incident SV wave is equal to the angle 
of reflection of the reflected SV wave, as shown in Fig. 2. 
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However, the angle of reflection of the reflected P wave is 
not equal to the angle of incidence of the incident SV wave. For 
a given value of b in Eq. (16), two values of the z-component of 
the slowness vector, S z , of equal magnitude but opposite sign 
(for a P wave travelling in the plane x - 0 in the transversely 
isotropic medium) can be obtained from the slowness surface for 
the P wave by substituting Eqs . (6) and (16) into Eq. (15). The 
positive z-component of the slowness vector corresponds to the 
reflected P wave and is denoted as S z ^^. Similarly, for a given 
value of b, there exists a positive z-component of the slowness 
vector, S z ( sv ), corresponding to the reflected SV wave. Accord- 
ingly, for any given value of b in Eq. (16), there exist a posi- 
tive z-component of the slowness vector for the reflected SV 
wave, S z < sv ), and a positive z-component of the slowness vector 
for the reflected P wave, S z ^). In fact, the z-components of 
the slowness vectors of the reflected SV and P waves, S z ^^^ anc * 
S z <*>, for an incident SV wave which determines the value of b in 
Eq. (16) and which travels in the plane x - 0 in the transversely 
isotropic continuum, can be obtained from the lengths of the per- 
pendicular lines between the abscissa representing the value of b 
and the intersections with the two sheets of the slowness sur- 
faces for the reflected SV and the reflected P waves, respec- 
tively, as shown in Fig. 3. It is apparent from Fig. 3 that the 
z-component of the slowness vector of the reflected SV wave, 
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S z ^ sv ), is greater than that of the reflected P wave, S z ( p ), for 
an incident SV wave travelling in the plane x - 0 in the unidi- 
rectional fiberglass epoxy composite [1] . Consequently, the 
angle of reflection of the reflected P wave 0p, defined similarly 
to Eq. (18) as 

0p - tan- 1 (S y ( p )/s z < p )>, (20) 

is larger than the angle of reflection of the reflected SV wave 
for an incident SV wave in the unidirectional fiberglass epoxy 
composite, from Eqs . (16), (18) and (20). 

The critical angle is defined as the angle of incidence of an 
incident SV wave at which the slowness vector of the reflected P 
wave becomes tangent to the plane boundary where the reflection 
occurs. For the specific case of an incident SV wave in the uni- 
directional fiberglass epoxy composite in Fig. 1, it has been 
shown above that the angle of reflection of the reflected P wave 
is greater than that of the reflected SV wave. Thus, the criti- 
cal angle phenomenon occurs for this composite. The critical 
angle for an incident SV wave in the unidirectional fiberglass 
epoxy composite shown in Fig. 1 is equal to 52°. The angle of 
reflection of the reflected P wave versus the angle of incidence 
of an incident SV wave in the unidirectional fiberglass epoxy 


17 



composite [1] is plotted from zero to the critical angle 52°, as 
shown in Fig. 4. 

4 . Amplitude Ratios of Reflected Waves to Incident Wave for 

Angle of Incidence Less Than Critical Angle 

It has been stated that an SV wave and a P wave will be 
reflected back into the semi -infinite transversely isotropic 
medium if the angle of incidence of the SV wave incident on the 
plane boundary is less than the critical angle. The stress 
boundary conditions for an incident SV wave whose slowness vec- 
tor is in the plane x - 0 can be rewritten, from Eqs . (4) and 
(7), as 


r yz 

r zz 


(I) 

(I) 


+ r < sv ) 
T r yz 


+ r (SV) 
r zz 


♦ 'yz' P) - 0 

+ 'zz (P) - 0 


( 21 ) 


where r and r zz^^ represent the shear and normal stresses 
associated with the incident SV wave; Ty Z (SV) and r zz (SV) 
represent the shear and normal stresses associated with the 
reflected SV wave; and ry Z ^ and r zz ^^ represent the shear and 
the normal stresses associated with the reflected P wave. 

The shear stress, Ty Z ^) , associated with an SV wave of unit 
amplitude travelling in the plane x - 0 and incident on the plane 
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boundary at the origin, as shown in Fig. 2, can be obtained from 
Eqs . (3) and (6) as 

r yz (I) - i W (C 55 S z ( I >P y < 1 ) + C 55 S y ( I )p z (D) (22) 

where P y (*) and P z ^ are the components of the unit vector of 
particle displacement of the incident SV wave along the y and z 
axes, respectively. Similarly, the normal stress, t zz ^) , 
associated with the incident SV wave of unit amplitude can be 
expressed, from Eqs. (3) and (6) as 

r zz (I > - i.(C 13 S y (I> Fy (I> + C 33 S Z <I>P Z <I>) (23) 


The shear stress, r yz ( sv ), and the normal stress, r zz; (SV) , 
associated with the reflected SV wave on the plane boundary at 
the origin are, from Eqs. (3) and (6) 

r yz (SV) . iw A(SV) ( c 55 S z (SV) Py (SV) + C 55 S y (SV)p z (SV) } (24) 

and 

r 2z <SV> - ia,A< SV >(C 13 s/ SV >P y < SV > + C 33 S 2 < SV >P z <SV)) (25) 
where p (SV) an(3 p (SV) are the components of the unit vector of 

y c 

particle displacement of the reflected SV wave along the y and z 
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axes, respectively; and A^V) is the amplitude of the reflected 
SV wave . 

The shear stress, Ty Z ^\ and the normal stress, r zz ^), 
associated with the reflected P wave on the plane boundary at the 
origin are, from Eqs . (3) and (6), 

r yz m . i^(P)(C 55 S z <P)p y (f) + C 55 S y (P)p 2 (P>) (26) 

and 

r zz < P > - i W A(P)(C 13 S y (P)p y <P) + C 33 S Z (P)P Z (P)) (27) 

where P y ^ and P Z ^ P ^ are the components of the unit vector of 
particle displacement of the reflected P wave along the y and z 
axes, respectively; and A^^ is the amplitude of the reflected P 
wave . 

Upon substitution of Eqs. (22) through (27) into Eq. (21), 
the stress boundary conditions for an incident SV wave of unit 
amplitude travelling in the plane x - 0 can be expressed as 
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( 28 ) 


S z ( I >Py< I > + S y (I)p z (I) + A(P)(S z (P)p y (P) + S y ( P >P z (P>) 

+ A (SV) (Sz (SV) Py (SV) + S y ( SV >P z ( SV )) - 0 

(C 13 s y ( 1 )p y (1) + c 33 s 2 (i)p z a)) + A (P)(C 13 S y (P>P y (P) 

+ C 33 S Z (P)P Z (P)) + A(SV)( Cl3 S y (SV) Py (SV) 

+ c 33 s z < sv )p z ( sv )) - 0 

The components of a unit vector of particle displacement 
along the y and z axes of an SV wave travelling in any plane 
containing the zonal axis of a transversely isotropic medium are 
given in [3], and when applied to the present case of the 
reflected SV wave travelling in the plane x - 0, can be expressed 
as 


P (SV) _ 


H SV ‘ ( c 33 * c 44) s z 


(SV)' 


H SV ’ ( c 33 - c 44> s z 


(SV)' 


( c 13 + c 44 )s<sv)s z (sv) 


1/2 


(29) 


P (SV) 
r z 


(C 13 + C 44 )S y (SV)s z (SV) 


H S V ' (^33 ’ C 44 )S z ( SV > i 


(C 13 + C 44 )S y (SV)s z (SV) 


1/2 
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where Hgy is defined as 




(Cn - C 44 )S(SV) 2 + (c 33 - C 44 )S y (SV) 2 


(Cn - C 44)Sy (SV) 2 + ( c 33 - C 44 )S Z (SV) 2 


- 4S y (SV) 2 s z (SV) 2 


( c ll - C 44 )(C 33 - C 44 ) 


- (C 13 + C 44 )‘ 


" 2 1 

1/2' 

► 

J J 



Similarly, the components of a unit vector of particle displace- 
ment along the y and z axes for the reflected P wave travelling 
in the plane x = 0 can be expressed as 


P (P) 


Hp - (C 33 - C 44 )S z ( p >‘ 


Hp - (C 33 - C 44 )S Z (P>- 


(Cl3 + C 44 )S y (P) Sz (P) 


1/2 


(30) 


P (P) 
r z 


(Ci 3 + c 44 )s y ( p )s z (p) 


( c l 3 + C 44 )s y ( p )s z ^ p ) 


Hp - ( c 33 ' C 44> s z 


(P)‘ 


1/2 
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I 

! 

I 

f 


where Hp is defined as 


Hp - |(C U - C A4 )S <P > 2 + (C 33 - C 44 )S z ( p ) 2 


+ < 


(Cll - C 44 )S ( p ) 2 + (C33 - c 44) s z ( p ) 2 


- 4S y ( p ) 2 S z ( p ) 2 


( c ll - C 44 )(C 33 - C 44 ) 


- (c 13 + c 44 )‘ 


' 2 1 

1/2' 

» 

J J 
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The amplitude ratios of the reflected SV and P waves, a( sv ) 
and A^ p ) , for an incident SV wave of unit amplitude whose angle 
of incidence is less than the critical angle can be determined 
from Eq. (28). For a given incident SV wave, the values of the 
components of the slowness vector along the y and z axes, S y ( p ) 
and S^ 1 ), and the values of the components of the unit vector of 
particle displacement along the y and z axes, P y '*' and P z ( p ), 
are defined as part of the specification of the incident SV wave. 
The values of the components of the slowness vector of the 
reflected SV wave along the y and z axes, S y (SV) anc i s z ^ V ), are 
determined from Eq. (16) and (17), respectively. The values of 
the components of the unit vector of particle displacement along 
the y and z axes for the reflected SV wave, P y ( s ^) and P Z ^V), 
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are obtained by substituting Sy( sv ) and into Eq. (29). 

The values of the coefficients for A^V) £q (28) are thus 
obtained. Since the y-component of the slowness vector of the 
reflected P wave, Sy( p ^ , is equal to that of the incident SV 
wave, Eq. (16), the value of the z-component of the reflected P 
wave, S z ( p ), is obtained by substituting Sy( p ) into the slowness 
surface for the P wave, Eq. (15). On substitution of the values 
of S z < p ) and Sy( p ) into Eq. (30), the values of the unit vector 
of particle displacement along the y and z axes for the reflected 
P wave, Py( p ) and P z ( p ), are determined. The values of the 
coefficients for A^ p ) in Eq. (28) are thus obtained. The values 
of the amplitudes of the reflected SV and P waves, and 

A( p ), are then obtained by solving Eq. (28) with the thus- 
determined values of the coefficients for A^ SV ) and A^ p ) in Eq. 
(28) for a given incident SV wave. By varying the angle of inci- 
dence of the incident wave. By varying the angle of incidence of 
the incident SV wave of unit amplitude and by repeating the 
procedures described above, the amplitude ratios of the reflected 
SV and P waves to the the incident SV wave are obtained as func- 
tions of the angle of incidence. For an incident SV wave in the 
unidirectional fiberglass epoxy composite as shown in Fig. 1, the 
amplitude ratios of the reflected SV and P waves to the incident 
SV wave versus the angle of incidence are plotted in Fig. 5 from 
zero to the critical angle 52°. 
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5 . Existence of Reflected Surface Wave for Angle of Incidence 


Equal to or Greater Than Critical Angle 

It is seen from Fig. 5 that at the critical angle 
52°, the amplitude ratio of the reflected SV wave to the incident 
SV wave is minus one and the amplitude ratio of the reflected P 
wave to the incident SV wave is -2.98 instead of being equal to 
zero. This suggests the possibility of the existence of a sur- 
face wave travelling parallel to the plane boundary for the angle 
of incidence equal to or greater than the critical angle. In the 
case of an incident SV wave in a semi- inf inite isotropic medium 
with an angle of incidence greater than the critical angle, there 
exist a reflected SV wave and a surface P wave whose amplitude 
decays exponentially with the perpendicular distance from the 
plane boundary [4] . 

Consider an SV wave travelling in the plane x - 0 in the 
semi- infinite transversely isotropic medium which is incident on 
the plane boundary, with the angle of incidence equal to or 
greater than the critical angle. Assume that there exists a sur- 
face wave travelling parallel to the plane boundary; that is, the 
x- component and the z- component of its slowness vector, S x and 
S z , vanish, whereas the y-component of its slowness vector, Sy, 
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does not vanish, as shown in Fig. 6(a). Let the reflected sur- 
face wave be represented as 


(u, v, w) - A< s )exp(-C 1 z)(P x < S ), P y < s >, P Z < S >) 

exp{ iw(Sy( s )y - t) } (31) 

where P x ^', P y (^) an( j p z (S) are the components of a unit vector 
of particle displacement of the reflected surface wave along the 
x, y and z axes, respectively; is the amplitude of particle 

displacement of the reflected surface wave; S y (^) is the y-com- 
ponent of the slowness vector of the reflected surface wave; and 
is an undetermined real constant which is equal to or greater 
than zero . 

Now, rearrange Eq. (31) into the following form: 

(u, v, w) - A( S >(P X < S >, P y < s >, P Z < S >) 

exp(iw[S y ( S )y + i(Ci/w)z] -t) (32) 

On substitution of Eq. (32) into Eq. (11), we find that 
P X (S) - 0 

[ c ll s y (S)2 - C 44 (C 1 / W )2 . p]P y (S) 
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+ (C 13 + C 44 )S y ( S )i(C 1 / w )P z ( S > - 0 


(33) 


(C13 + C 44 )S y ( S )i( Cl / W )P y (S) 

2 

+ [C 44 S y (S) - C 33 <C 1 /w ) 2 - p]P z ( s ) - 0 

The condition for the existence of non- trivial solutions for 
P y (^) and P z ^) is expressed by setting the determinant of the 
matrix of the coefficients of P y (S) anc j p z ( s ) in Eq. (33) equal 
to zero: 

2 

[CuS y (S) - ^(Ci/w ) 2 - p] 

( C 13 + C 44 )S y ( S )i(C 1 /u)) 

(34) 


(C 13 + C 44 )S y ( S >i(C 1 / W ) 


[c 44 s y (S) - C 33 (C 1 / W ) 2 - p) 


Expanding Eq. (34), wt obtain 


c 33 c 44 


((Vo)' 


C 33 (p - CllSy^ ) 


(35) 


+ C 44 (p - C 44 S y ^ S ^ ) + (Cjj + C 44 ) 2 s y ^ s ^ 


(C]_/w ) 2 


2 2 

+ (p - c n s y < s > ){p - c 44 s y < s ) ) - 0 


Eq. (35) is a quadratic equation in (C^/w) 2 ; designate its 
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coefficients as follows: 


A “ c 33 c 44 

2 2 

B - C 33 (p - C n S y ( S > ) + C 44 (p - C 44 S y ( S ) ) 

+ (C 13 + C 44 ) 2 S y (S> 2 (36) 

c - (p - c n s y < s > 2 )(p - c 44 s y < s ) 2 ) 

It is apparent from Eq. (36) that A is always positive. 

By definition, at the critical angle, the angle of reflection 
of the reflected P wave 0p is equal to 90° [4] ; accordingly, 

S z <*> in Eq. (20) is equal to zero. Therefore, the value of b in 
Eq. (16) corresponding to the critical angle is thus determined 
by setting S x and S z in Eq. (15) equal to zero, and is found to 
be (p/ C]_]_)^/ 2 . As stated following 2q. (19), the angle of 
reflection of the reflected SV wave 8§y is equal to 90° when the 
SV wave is tangentially incident on the plane boundary; accord- 
ingly, s z < sv > in Eq. (18) is equal to zero. Therefore, the value 
of b in Eq. (16) corresponding to the incident SV wave whose 
slowness vector is parallel to the plane boundary and is in the 
plane x - 0 is determined by setting S x and S z in Eq. (14) equal 
to zero, and is found to be (p/C 44 )^/ 2 . The necessary condition 
for the occurrence of the critical angle phenomenon for an SV 
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wave travelling in the plane x - 0 is that the value of 
(P/C 44) 1 / 2 has to be greater than that of (p/Cn)l/2 [ 5 ] Take 
numerical values of (p/C 44) ^/ 2 and (p/Cn )^ 2 for 
unidirectional fiberglass epoxy composite shown in Fig. 1 from 
[ 1 ] as follows: (p/ C j_j_)V 2 _ 0.418(s/km) and (p/ 044 )^/^ - 
0.647(s/km), as shown in Fig. 3. Therefore, these values satis- 
fying the necessary condition on the occurrence of the critical 
angle phenomenon for an incident SV wave in the unidirectional 
fiberglass epoxy composite. 

It has been shown [5] that when the angle of incidence of an 
incident SV wave travelling in the plane x - 0 is equal to or 
greater than the critical angle, the value of b in Eq. (16) lies 
between (p/Cj_j _)^/ 2 and (P/C 44 ) , that is; (p/C;q )^/ 2 < b < 
(P/C 44 ) 1 / 2 , as shown in Fig. 3. In addition, it follows from Eq. 
(5) that the y-component of the slowness vector of the reflected 
surface wave, Sy(^) , is equal to b in Eq. (16). As a result, the 
value C in Eq. (36) is equal to or less than zero when the angle 
of incidence of an incident SV wave travelling in the plane x - 0 
is equal to or greater than the critical angle; that is, C < 0. 
Accordingly, the value of B 2 - 4AC defined in Eq. (36) is equal 
to or greater than zero for the angle of incidence equal to or 
greater than the critical angle; that is, 
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B 2 - 4AC - B 2 + 4A|C| > 0 


(37) 


Thus, the possible solutions to the quadratic equation in (q/w) 2 
in Eq. (35) .are 


0 -B ± (B 2 - 4AC) 1 / 2 

( Cl / W ) 2 — (38) 

2A 

However, it follows from Eq. (37) that the value of B 2 -4AC is 
equal to or greater than that of B 2 for any angle of incidence 
equal to or greater than the critical angle. Consequently, the 
solution to (q/a>) 2 , from Eq. (38), is 

0 -B + (B 2 - 4AC) 1 /2 

(Cj/w) 2 0 (39) 

2A 

Accordingly, there exists a real constant q given in Eq. (31) 
when the angle of incidence of the incident SV wave travelling in 
the plane x - 0 is equal to or greater than the critical angle, 
and is expressed from Eq. (39), as 


-B + (B 2 - 4AC) 1 / 2 


2A 


> 0 


(40) 


The existence of a real constant q in Eq. (40) confirms the 
existence of the surface wave assumed in Eq. (31). 


It is therefore concluded that when the angle of incidence of 
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an incident SV wave travelling in the unidirectional fiberglass 
epoxy composite as shown in Fig. 1 is equal to or greater than 
the critical angle, only an SV wave is reflected back into the 
medium, and the reflected P wave degenerates into a surface wave 
which travels parallel to the plane boundary and whose amplitude 
decreases exponentially with perpendicular distance form the 
plane boundary, as shown in Fig. 6(b). 

6 . Balance in Energy Flux Normal to Plane Boundary 

The balance in energy flux normal to the plane boundary must 
be satisfied [2] . For the angle of incidence less than the crit- 
ical angle, the balance in energy flux normal to the plane bound- 
ary z “ 0, as shown in Fig. 2, is expressed as [2] 

F Z (I) + F z ( p ) + F Z < SV > - 0 (41) 

where F z ^ p \ F Z ^V) and F z ^ p ) are the z -components of the energy 
fluxes of the incident SV wave, the reflected SV wave and the 
reflected P wave, respectively. 

The z- component of the energy flux of an incident SV wave of 
unit amplitude travelling in the plane x - 0 is [6] 
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) 



( 42 ) 


F z (D - u 2 (C44 p y (I> 2 s 2 (I) + C 13 P y U)p 2 (I)s y U) 

+ c 44 p y < I >p z < I )s y a) + c 33 p 2 (i) 2 s z a)) 

Similarly, the z-components of the energy fluxes of the reflected 
SV and the reflected P waves are [6] 

P 2 (P) -a(P) 2 „2 ( c 44 p (P) 2 s 2 (P) + c 13 P v (P)p z (P)s y (P) 

(43) 

+ G 44 P y <P)p 2 <P)s y <P) + C 33 P 2 (P) 2 S 2 (P)) 

and 

Fz (SV) » A ( SV ) 2 W 2 (C 44 P y (SV) 2 s z (SV) + C 13 P y ( SV )p z ( SV ) Sy ( SV ) 

+C 44 P y (SV)p z (SV) Sy (SV) + c 33 p z ( sv > 2 s 2 ( sv )) (44) 

Since the frequency term o> 2 is common to Eqs. (42), (43) and 
(44), the balance in energy flux normal to the plane boundary, 

Eq. (41), is not affected by assuming the value of the frequency 
u to be equal to unity. Accordingly, subsequent calculations of 
the values of the z-components of the energy fluxes of the inci- 
dent SV wave and the reflected SV and P waves, F z ^), F Z ^V) an<3 
F z ^), are done by assuming the radian frequency w to be equal to 
one . 
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The value of the z-coraponents of the energy fluxes in Eqs . 
(42) ,(43) and (44) are obtained similarly to the calculations of 
the amplitude ratios of the reflected P and SV waves to the inci- 
dent SV wave for the angle of incidence less than the critical 
angle. For a given incident SV wave of unit amplitude, the 
y-components of the slowness vector and of unit vectors of par- 
ticle displacement, , Sy( SV ), Sy( p ) , Py^), Py^ SV ) and 

Py^), and the z-component of the slowness vectors and of the 
unit vectors of particle displacement, S z ^), S z ( sv ), S z ^^, 
P z ^), P z ^^ and P z ^) of the incident SV wave and the reflected 
SV and P waves are determined first. Then combining the ampli- 
tudes of the reflected SV and P waves, and , with the 

values of Sy (D , Sy (SV) , jy^, S Z <D, S Z <SV), s z <*>, P y <I>, 
Py(SV), Py(P), P z ^^, P z ^v) and P z ^' , the values of the 
z-components of the energy fluxes of the incident SV wave 
and the reflected SV and P waves, , F z ^ v ) and F z ^\ are 

thus obtained from Eqs. (42), (43) and (44). 


For an SV wave of unit amplitude travelling in the plane x - 

0 in the unidirectional fiberglass epoxy composite shown in Fig. 

1 incident on a plane boundary, the z-components of the energy 
fluxes of the reflected SV and P waves and the energy flux of the 
incident SV wave are plotted in Fig. 7 from zero to the critical 
angle 52°, with the value of the frequency w in Eqs. (42), (43) 
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and (44) equal to one. The balance in energy flux normal to the 
plane boundary, Eq. (41), is also checked and is shown in Fig. 7. 

When the angle of incidence of an incident SV is equal to or 
greater than the critical angle, only an SV wave is reflected 
back into the medium with its amplitude ratio to the incident SV 
equal to minus one, and the reflected P wave degenerates into a 
surface wave whose energy flux is parallel to the plane boundary 
[6], As a result, the energy flux normal to the plane boundary 
due to the surface wave is equal to zero. The balance in energy 
flux normal to the plane boundary z - 0, as shown in Fig. 6(a), 
for the angle of incidence equal to or greater than the critical 
angle is expressed as [2] 

F (I) + F (SV) _ o (45) 


The relationships between the components of the unit vector 
of particle displacement of an incident SV wave travelling in the 
plane x - 0 along the y and z axes, Py(^) and P z ^^, and those of 
the reflected SV wave can be obtained by substituting Eqs . (16) 
and (17) into Eq. (29) and can be expressed as 


P (I) 


P (SV) 

r y 


p 2 d) - -p z (sv) 


(46) 
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On substitution of Eqs . (16), (17) and (46) into Eqs . (42) 
and (44) , it is found that 

f z CD - -F z < sv > (47) 

Accordingly, the balance in energy flux normal to the plane bou- 
ndary, Eq (45), is satisfied when the angle of incidence of an 
incident SV wave is equal to or greater than the critical angle. 
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ACOUSTO-ULTRASONIC NON-CONTACT INPUT-OUTPUT CHARACTERIZATION OF 


UNIDIRECTIONAL FIBERGLASS EPOXY COMPOSITE PLATE 

It has been shown that the unidirectional fiberglass epoxy 
composite shown in Fig. 1 may be modelled as a homogeneous trans- 
versely isotropic continuum. For the axes shown in Fig. 1, the 
isotropic plane of the equivalent continuum lies in the midplane 
of the plate [1] . A cartesian coordinate system (x, y, z) is 
chosen so that the x-y plane is the isotropic plane; thus, the 
upper and lower faces are at z - h/2 and z - -h/2, respectively, 
where h is the plate thickness. The properties of the equivalent 
continuum model of the unidirectional fiberglass epoxy composite 
plate to be considered are [1] 

h - 0.1 m 

C n - 10.581 x 10 9 N/m 2 

C 13 - 4.67 x 10 9 N/m 2 (48) 

C 33 - 40.741 x 10 9 N/m 2 
C 44 - 4.422 x 10 9 N/m 2 
C 66 - 3.243 x 10 9 N/m 2 
p - 1850 kg/m^ 

Non- contact transmitting and receiving transducers are 
located on the same face of a fiberglass epoxy composite plate 
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specimen, as shown in Fig. 8. The unidirectional fiberglass 
epoxy composite plate specimen is considered as a plate of thick- 
ness h and of infinite planar (x-y) extent. The input electrical 
voltage to the transmitting transducer is V^(t) and the output 
electrical voltage from the receiving transducer is V Q (t) where t 
represents time. The transmitting transducer converts an input 
electrical voltage into a stress, whereas the receiving trans- 
ducer converts a displacement associated with stress waves trav- 
elling in the plate into an output electrical voltage. In the 
following analysis, only the SV waves are traced. The SV waves 
which are generated by the transmitting transducer located above 
point 0 experience multiple reflection at each face of the plate, 
and then reach the receiving transducer located above point M, as 
shown in Fig. 9. Since the isotropic plane lies in the midplane 
and is parallel to both the top and bottom faces where the mul- 
tiple reflections occur, the angle of reflection of the reflected 
SV wave is equal to the angle of incidence of an incident SV wave 
for each reflection at each face of the plate. Accordingly, the 
SV waves travelling from the input 0 to the output M may be con- 
sidered as waves propagating in a semi - inf inite transversely iso- 
tropic medium and travelling to point M' as if there were no 
bottom face, as shown in Fig. 9. 
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1 . Delay Time and Phase Velocity 


Let the input 0 and the output M lie in the y-z plane. 

Assume the number of reflections at the bottom face experienced by 
the SV wave in travelling from the input 0 to the output M is n, 
as shown in Fig. 9. With respect to the z axis, the angle of 
incidence of the SV wave at each face of the plate is 8 , and the 
total distance travelled by the wave is R^. From the geometry in 
Fig. 9, 


8 - tan _1 (i/2nh) (49) 

where 2 is the separation distance between the input 0 and the 
output M, 

Rn = 2 / sin# (50) 

The delay time t n for the SV wave to reach the receiving trans- 
ducer is 

t n “ Rn/Cl(*) (51) 

where C^(0) is the directionally -dependent phase velocity of the 
SV wave. The phase velocity of an SV wave in the unidirec- 
tional fiberglass epoxy composite plate is [1] 

C].(0) - [(C 44 + C n sin 2 0 + C 33 cos 2 0 - yl)/2p] 1/2 (52) 
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where 


2 

E ” - C 44 )sin ^0 + (C 44 - 

+ 4(C^3 + C44)^sin^Scos^fl ; 

c ll> c 13< C 33> c 44 and P are given by Eq. (48). 

The phase velocity as a function of the angle of incidence 8 
is shown in Fig. 10. 

The delay time is then computed when the number of reflec- 
tions n at the bottom face of the plate is equal to 10, 100, 300 
or 500. The numerical results are shown in Fig. 11 where the 
delay time t n is plotted as the ordinate, and the dimensionless 
separation i/h is plotted as the abscissa, for values of zero to 
300. 

2 . Displacements Detected bv Receiving Transducer 

The displacements detected by the non-contact receiving 
transducer above point M, radiated by the non-contact transmit- 
ting transducer, are assumed to be equivalent to the displacement 
at point M' associated with the SV wave propagating in a semi- 
infinite transversely isotropic medium. Except for the reflection 
coefficients at each face (to be discussed later), the displace- 
ment is computed as if there were no bottom boundary, as shown in 
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Fig. 9. The displacement at point M' is approximated by the 
far- field asymptotic solution for large of an infinite trans- 
versely isotropic medium subjected to a harmonic point force. 

Consider an infinite transversely isotropic medium in which 
the z axis of a rectangular cartesian coordinate system 0(x, y, 
z) is the zonal axis of the medium, and the x-y plane coincides 
with isotropic plane, as shown in Fig. 12. The equations of 
motion including the body force are [7] 

r xx, x + r xy,y + r xz,z + PX ” pu.tt (53) 

r xy, x + r yy,y + r yz , z + pY — pv,tt (54) 

r xz, x + r yz,y + r zz,z + pZ “ pw.tt (55) 

where r rs (r, s - x, y and z) are the normal (r - s) and shear 
(r * s) stresses with respect to the chosen coordinate system 
0(x, y, z) ; u, v and w are displacement components of a point in 
the medium along the x, y and z axes, respectively; X, Y and Z 
are the components of the body force along the x, y and z axes, 
respectively; p is the density; t is time ; and denotes 
partial differentiation with respect to the variable which fol- 
lows . 
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Combining Eqs . (1) and (55) gives 


C]_3 + C44 C44 Caa 

r.tt - (— ) A,Z2 + (- 2 =-) r,xx + (-==-) r.yy 

P p p 

C33 

+ ( — — ) r.zz + z, z 

p 


( 56 ) 


where T and A are given by [6] 


T - w,z and A - u,x + v,y . 

By differentiating Eq . (53) with respect to x and Eq. (54) with 
respect to y, we find, on addition of the resulting equations and 
using the appropriate stress-strain relation, Eq . (1), that 


A t-f 


■ (r,xx + r,yy) + A, 22 


+ 


C 11 

p 


(A, xx + A,yy) + X,x + Y,y 


(57) 


For a harmonic point force at the origin, the body forces may 
be taken of the form [7] 


X - X Q S (k)S (y) <S (z )e~ iwt 

Y - Y Q S (x) 8 (y) 6 (z)e‘ iwt (58) 

Z - Z 0 6(x)6(y)6(z)e- iwt 


where S( r) (r - x, y and z) is the Dirac delta function. 
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Express r and A as threefold Fourier Integrals [7] as 


follows : 


r(x,y,z,t) 


f ( Sjj , Sy , Sj 


, t) 


exp{iw(S x x + Syy + S z z - t))dS x dSydS z 


(59) 


A(x,y,z,t) 


A(S X ,Sy,S Z ,t) 


exp(iw(S x x + Syy + S z z - t))dS x dSydS z 


where 


r(S x ,S y ,Sz,t) 


1 / 8 *' 


r(x,y,z,t) 


exp{ico(-S x x - Syy - S z z + t)}dxdydz. 


A(S X » Sy , S. 


; ,t) - 1/8*' 


A(x,y , z , t) 


exp{iw(-S x x - Syy - S z z + t))dxdydz. 
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Similarly, X,x, Y,y and Z,z can be expressed as threefold Fourier 
integrals 


X,x 


iwS, 


X exp{iw(S x x + Syy + S 2 z - t)}dS x dSy 


dS, 


(60) 


00 


Y,y = iwSy 


Y exp{iw(S x x + Syy + S z z - t))dS x dSydS z 


(61) 


00 


Z,z — iwS z 


Z exp{iw(S x x + Syy + S z z - t))dS x dSydS z 


(62) 


where 


X 


1/8jt 3 


X 0 5(x)6(y)5(z)e- iwt 


exp{iw(-S x x - Syy - S z z + t)dS x dSydS z - X 0 /8tt 3 
and, similarly, Y - Y Q /8ir 3 ; Z - Z 0 /8 tt 3 . 
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Substitution of Eqs . (58) through (62) into Eqs . (56) and (57) gives 


A = 


C 44 /p(S x 2 + Sy 2 ) + C 33 /pS z 2 -1 
8?r 3 a>H(S x , Sy, S z ) 


(iS x X Q + iS y "i 


- (Cj3 + C 4 4 )/P ( s x 2 + s y 2 ) 

8ff 3 wH(S x , Sy , S z ) 1 z 0 


- (0^3 + C44 ) / p S z 

— i 3 ^ — (is x 0 + is v Y 0 ) 

8ff 3 wH(S x ,S y ,S z ) y 0 


C n /p(S x 2 + Sy 2 ) + C 33 /pS z 2 -1 

T J iS z Z-, 

87r 3 o)H(S x ,Sy,S z ) 


where 


H(S x ,S y ,S z ) 


^ S 2 2 + 

IP P 


(S x 2 + Sy 2 ) - 1 


2 , c 2\ , _^33c 2 


r (S* + S v ^) + 


V -1 


( C — ■ 2 S Z 2 (S X 2 + S 2 ) 
P 


In fact, H(S X , Sy, S z ) - 0 represents the two sheets 
surface, one for a P wave and one for an SV wave [8] 
result of Eqs. (63) and (64), 


o) 

(63) 


(64) 


(65.) 


of slowness 
As a 
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r(x,y,z,t) 


-(C13 + C 44 )/p S z 2 ] 


[ 8jt j wH(S x , Sy ,S z ) 


(iS x X 0 + iSyY Q ) 


+ f c H /p(s x 2 + s y 2 > + c 33/^ s z 2 - 1 

8jr^wH(S x ,Sy, S z ) 


exp{iw(S x x + S y y + S z z - t) } 


i- s z z o 


( 66 ) 


A(x , y , z , t) 


C 44 /p(S x 2 + s y 2 ) + c 33 a>s z 2 -1 ' 

8AH(S x ,S yl S z ) 


(iS x X G + iS y Y 0 ) 


(Cj 3 + C 44 )//? (S X 2 + Sy 2 ) 
— 8 jt 3 wH(S x , S y ,S z ) ^ 


iS z Z 0 


( 67 ) 


exp(iw(S x x + S y y + S z z - t) } 


The asymptotic solution at a large distance from the point 
force is obtained by applying the theory of residues, the method 
of stationary phase, and the radiation condition [8] as 


r(x,y,z,t) 


w N 

— — zy n 

2 jtR n -1 


-(Ci 3 +C 44 ) 2 

S zn (^ s xn^o + ^^yn^o^ 


C 11 2 2v f33 s 2 . x 

( s xn + b yn ' + b zn 1 
l P J p 


^zn 2 o 


( 68 ) 


exp{ iw(S xn x + Sy^y + S z ^z - t) ) 
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A(x,y,z,t) 


u> 

2?rR 


n rrc 44 


^xn + ^yn ) + ' 


'33 


S - 1 
b zn i 


^ S xn X o +: *- S yn^o^ + 


( c 13 +c 44 ) 


( s xn 2+ s yn 2) iS ™ Z 


°} 


exp{iw(S xn x + Sy n y + S zn z ■ t) } 


( 69 ) 


where R is the distance from the origin 0 where the point forces 
are applied to the location of interest Q in the medium, as shown 
in Fig. 12; (S xn , Sy n , S zn ) are points on the slowness surface 
for an SV wave where the normal is parallel to the OQ direction; 

N is the total number of points (S xn , Sy n , S zn ) are points over 
which the summation must be performed ; A n is the amplitude 
coefficient and is given by 


_ f H,S X 2 + H,S y 2 + H,S Z 2 

An l Kl J 


1/2 


(70) 


where Kjj is given as 


kn = 2 [H,S z 2 (H,S x S x H.SySy - H,S x S y 2 ) 

+ 2H,S X H,S y (H,S x S z H,S y S z - H,S x S y H,S Z S Z )] 


where 2 denotes the sum with respect to cyclic permutation of S x , 
Sy and S z ; the symbol | | denotes "the magnitude of"; A n is 

evaluated at points (S xn , Sy n , S zn ) on the slowness surface for 
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an SV wave where the normal is parallel to the OQ direction, 
exclusive of those singular points (S xn , S yn , S zn ) resulting in 
K n in Eq. (71) equal to zero; and A n is the phase constant and is 
determined as follows: A n - 1 if K n > 0 or - i if K n < 0. 


The displacement components along the x, y and z axes, u, v 
and w due to an SV wave can be obtained by direct integration of 
the definition of T and A in Eq . (56), and are given as [8] 


u 


iw(S v z + S v z ) 


iu(S x 2 + S y 2 ) i 


(72) 


1 

w = r 

iwS z 

Substitution of Eqs . (68) and (69) into Eq. (72) gives the 
asymptotic solutions at a large distance to the displacement com- 
ponents along the x, y and z axes, u, v and w, as follows: 


-i N 

u ~ — SV^n 

2?rR n~l 


J xn 


l (Sxn 2 + Syn^ ) L 


c 44 
l P 


( s xn 2 + s yn 2 ) + ' 


J 33 


'zn 


- 1 


(iS xn 2 X 0 + iSy n 2 Y 0 ) + ~ (Cl 3 +C - 4 ) (S xn 2 + S yn 2 ) iS zn 2 Z 


exp(iw(S xn x + S yn y + S zn z - t) ) 


(73) 
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V ~ 


-i N 


2 A n X n T 

2nR n-1 ^ n (S xn 2 + S^ 2 ) 


44 (S 2 +S 2 )+ S 2 - 1 

^xn a yn > + s zn • L 


l P 


^ s xn X o + ^yn + 


- ( c 13 +c 44) 


+ S yn^) ^zn Z o 


exp(iw(S xn x + S^y + S zn z - t) } 


(74) 


1 N 

w ~ — XA n A n 
2 ttR n=l 


-(C 13 +C 44 ) 


s zn^ iS xn x o + ^ s yn Y o^ 


+ < — ( s xn 2 + S yn 2 ) + ~ S zn 2 * DZ e 


P ' P 

exp(iw(S xn x + S^y + S zn z - t) ) 


(75) 


3 . Directivity Functions 


The shear stress Ty Z and the normal stress r zz associated 
with the SV waves reaching the point M' in Fig. 9 are used to 
study their associated directivity functions. The asymptotic 
shear stress Ty Z and the asymptotic normal stress r zz are 'given 
py substituting Eqs. ( 68 ) through (72) into Eq. (1) as 


zz 


w n r 

~ — XA n A n j c 33 
2?rR n=l ^ 


-(C 13 +C 44 ) 9 

s zn ( iS xn x o +: '-^yn Y o^ 
P 


_ii_(s 2 + g 2 \ [ 33 g 2 _ ^ 
^xn + a yn > + s zn 1 

P J p 


^zn^o 
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'13 


c 44 


^xn + s yn ) + 


c 33 


l P 

-(c 13 +c 44 ) 


zn 


¥ 


1 ( IS^Xq+IS^Yq) 


^xn + ^yn ) ^zn^o 


1 } 


exp { iw ( S xn x + S^y + S zn z - t) 


(76) 


yz ~ 


C 55 u> 

2nR 


N f r-(Ci3 + C 4 4) - 

^^n |V[ ' s zn (iS xn X 0+ i Syn Y 0 ) 


(_£ii(S 2 + g 2\ + ^33 

^ ( ' a xn ^ a yn > + “zn -/--“o 

P P 


S zn Z - l)iZ f 


S^S 




zn 


o Z . O z 

a xn + a yn 


r ^44 

p 


t . ^xn^ + Syn^) + ^ zn Z ' (^xnX 0 + i.Sy n Y 0 ) 


-(C13+C44) 9 9 

+ < S xn + s yn 2) iS zn z o 


exp(iw(S xn x + S^y + S zn z - t) ) 


(77) 


The directivity functions associated with the normal stress 
r zz in Eq. (76) and the shear stress Ty Z in Eq. (77) will be 
evaluated for the case of S x - 0. Due to the axial symmetry with 
respect to the zonal axis, the z axis, of the transversely iso- 
tropic medium, the values of the directivity functions thus 
obtained hold for all values of S x . 


Consider the case of the point force acting along the y 
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direction only; that is, Y 0 * 0 but X 0 - Z Q - 0 in Eq. (58). 

Y 

The directivity function Dy Z of the shear stress r y Z 
associated with an SV wave is evaluated from Eq. (77) along the 
arc y^ + - 1 in the positive y-z quadrant by setting Y Q - R = 

l; 



tjC 55 

2n 


= Wzn( 
n-1 





1 ) 


exp(iw(S yn y + ^zn z ) ^ 


(78) 


where the summation 2 is done vectorially in a complex plane over 

each type of wave passing a given point (0, y, z) on the arc y^ + 

- 1; each wave corresponds to a point (0, S^, S zn ) on the 

slowness surface of an SV wave at which the normal is parallel to 

the radius vector of the given point (0, y, z). Similarly, the 

y 

directivity function D zz of the normal stress r zz is, from Eq. 
(76), for Y 0 - R - 1; 




N 

^n^n^yn 


c 13 c 44 „ 2 

p 3 ™ 


c 33 c 44 „ 2 

b zn 
P 



exp{iw(S yn y + S zn z) ) 


(79) 


Next, consider the case of the point force acting along the z 
direction only; that is, Z Q * 0 but X Q - Y Q - 0 in Eq. (58). 


Z 
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The directivity function Dy Z of the shear stress Ty Z 

associated with an SV wave is evaluated from Eq. (77) along the 
9 9 

arc y^ + z z - 1 in the positive y-z quadrant by setting Z 0 - R = 

l; 



^55 

2?r 


N 

Wyn 
n-1 J 


C 11 
. P 


„ 9 C 33 - C 44 - c 13 

N + " 



expfiwCS^y + S zn z)} 


( 80 ) 


Similarly, the directivity function D zz of the normal stress r zz 
is, from Eq. (76), for Z Q - R - 1; 



N 

^ ^n^n^zn 
n-1 


c ll c 33 ~ c 13 c 44 ~ c 13 2 
P 



+ 



P 




exp{iw(S yn y + S zn z)} 


(81) 


Finally, consider the case of the point force acting along the 
x direction only; that is, X Q + 0 but Y 0 - Z 0 - 0 in Eq. (58). It 
follows from Eqs. (76) and (77) that the values of the shear 
stress Ty Z and the normal stress r zz associated with an SV wave 
travelling in the plane x - 0 are equal to zero. 


The polar diagrams for the directivity 

Y 

stress Ty Z and the normal stress r zz (Dy Z , 


functions of the shear 
D Iz> D yz and D zz given 
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in Eqs . (78) through (81) ) associated with the SV waves propa- 


gating in the unidirectional fiberglass epoxy composite as shown 
in Fig. 1 are obtained by substituting Eq. (48) into Eqs. (78) 
through (81) at frequencies of 0.75, 1.50 and 2.25 MHz. 

Numerical results are shown in Figs. 13 through 24 where the 
angle of incidence 6 given in Eq. (49) is used to determine the 
direction in which the values of points (0, Sy n , S zn ) on the 
slowness surface, where the normal is parallel to the given 
direction, is thus obtained. 

4 • Assumptions of Transducers 

The non- contact transmitting transducer in Fig. 8 is assumed 
to transform an electrical voltage into a uniform stress; how- 
ever, the non-contact receiving transducer in Fig. 8 transforms a 
displacement into an electrical voltage. The approach below is 
similar to that given in [9]. Referring to Fig. 8, if an input 
voltage of amplitude of V and frequency u> is applied according to 

V £ (t) - Ve' i<Jt . (82) 

The stress a that is introduced into the specimen plate by the 
non-contact transmitting transducer is 
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cr(t) - F 1 (w)Ve' i < wt + ^l) 


(83) 


where F^(w) is the transduction ratio for the non-contact 
transmitting transducer in transforming a voltage to a stress 
and <f>^ i- s a phase angle. In Eqs . (82) and (83), the complex 
harmonic character of the signals is expressed in the complex 
notation where i - 7-1 and only the real part of these and 
subsequent equations should be considered. Thus, the amplitude T 
of the applied load is defined as 

T - F 1 (w)V (84) 

Similarly, if a stress wave producing a displacement d of 
amplitude D and frequency w that is detected by the non-contact 
receiving transducer is defined as 

d(t) - De‘ i(Jt (85) 

the output voltage from the non-contact receiving transducer, see 
Fig. 8, is 

V Q (t) - F 2 (w)De* i ( wt + W (86) 

where F 2 (w) is the transduction ratio for the non-contact 
receiving transducer in transforming a displacement to a voltage, 
and 4>2 is a phase angle. Thus, the amplitude V' of the output 
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electrical voltage is 


V' - F 2 0)D (87) 

The characteristics of Fi(w)F 2 (w) are unknown except that the 
dimensions of the product F^(w) are [kg/m^ • sec^ ] . 

5 . Steadv-State Output Voltage Amplitude due to Multiple Wave 
Reflections in Plate 


Since the SV wave traced in the unidirectional fiberglass 
epoxy composite plate specimen shown in Fig. 9 is travelling in 
the y-z plane, it follows from Eqs . (73), (74) and (75) that only 
the displacement components along the y and z axes, v and w, are 
detectable at the point M' . Consider first the point force 
acting along the y direction only; that is Y 0 * 0 but X Q - Z Q - 0 
in Eq. (58) . 

y 

The amplitude of the y-component displacement D v evaluated at 
the point M' can be obtained from Eqs. (50) and (74) as 


^l^ s yn> s zn) ,, 
u v „ *o 


where 


•^l^ s yn’ s zn^ 


Rn 


1/2tt 


( 88 ) 


N 


2 • A n*nt< c 44A’> s yn 2 + ( c 33^) s zn 2 
n-1 J 


- 1] 


exp { iw ( Sy n y + S 2n z)} 
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and is evaluated along the arc y 2 + x 2 - 1. Similarly, the 

Y 

amplitude of the z-component displacement D w evaluated at the 
point M' can be obtained from Eqs. (50) and (75) as 


Y f 2( s yn> s zn> 


w 


yn’ 

Rn 


( 89 ) 


where 

f 2 ( S yn’ S zn) ’ V2* 


N 


2 + c 44)/^l s yn s zn 

n-1 


exp{iw(S yn y + S 2n z)} 


Next, consider the point force acting along the z direction 
only; that is Z Q r* 0 but X 0 - Y 0 - 0 in Eq. (58). 

2 

The amplitude of the y- component displacement D v evaluated at 
the point M' can be obtained from Eqs. (50) and (74) as 


f 2^ S yn> s zn> „ 
Dv ^ Z 0 


(90) 


Similarly, the amplitude of the z-component displacement D w 
evaluated at the point M' can be obtained from Eqs. (50) and (75) 
as 


„Z f 3< s yn> s zn> „ 
«n ° 

where 


(91) 


■^3^yn’ ®zn^ 


1/2jt 


S A n A n [(C 11 /p)S 2 + (C 33 /p)S zn 2 - 1] 
n-1 
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exp { iw(Sy n y + S zn z)} . 

Finally, consider the point force acting along the x direc- 
tion only; that is, X Q * 0 but Y Q - Z Q - 0 in Eq. (58). It 
follows from Eqs . (74) and (75) that the amplitude of the y-com- 
ponent displacement and the amplitude of the z- component dis- 
placement vanish at the point M' . This is due to the fact that 
the SV wave is travelling in the y-z plane. 


According to Eqs. (88) through (91), the amplitude of a dis- 
placement component at the point M' , denoted as D^' , can be 
expressed in the following form as 




fi(S 


ZE1 


S zn^ 


*n 


(92) 


where T is the applied point force and is equivalent to either Y Q 
or Z 0 in Eqs. (88) through (91); and S zn ) (i - 1, 2 or 3) 

is determined from one of Eqs. (88) through (91), depending on 
which displacement component is measured and along which direc- 
tion the point load is applied. 


However, with the bottom boundary present, the wave is 
reflected a total of (2n -1) times, as shown in Fig. 9. Thus, 
the amplitude of displacement at the point M is and is 
expressed as 
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(93) 


2n- 1 


% “ Qsvsv d M' 


where Qgysv i- s t ^ le amplitude ratio of the reflected SV wave to 
the incident SV wave, as shown in Fig. 5, and is a function of 
the angle of incidence 8 defined in Eq. (49). Here, we have 
ignored the effects of mode conversion. 


The amplitude of the output voltage form the non-contact 
receiving transducer is V' and can be obtained by substituting 
Eqs. (92) and (93) into Eq. (87) as 

2n-l 

v , _ F 2 (w)Q svsv f i (S yn , S zn ) ^ 

Rn 

Substitution of Eq. (84) into Eq. (94) gives 


2n- 1 

F 1 (u>)F2(^)QsvsV f j( S yn' S zn^ v 

R n 


(95) 


Introducing the SV wave attenuation constant a of the unidirec- 
tional fiberglass epoxy composite and a possible electrical sig- 
nal amplification factor K, Eq. (95) can be written as 

2n- 1 e' oRn 

V' - KF 1 ( W )F 2 ( W )Q SVSV fi(S S zn )V (96) 

Rn 

Eq. (96) gives the output voltage amplitude from the non-contact 
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receiving transducer due to an input voltage amplitude V at the 
non- contact transmitting transducer when the SV wave path has 
undergone n reflections from the bottom face of the unidirec- 
tional fiberglass epoxy composite plate specimen, as shown in 
Fig. 9. 
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DISCUSSION AND CONCLUSION 


In the acousto-ultrasonic input-output characterization of the 
unidirectional fiberglass epoxy composite plate, the angle of 
reflection of the reflected SV wave is equal to the angle of 
incidence of the incident SV wave for each reflection at either 
the top or the bottom face of the plate. This is due to the fact 
that the isotropic plane is parallel to both faces of the plate. 
However, if the parallelism between the isotropic plane and the 
plane boundaries where reflection occurs does not exist, the 
angle of reflection is not equal to the angle of incidence. In 
such a case, the use of a semi -infinite transversely isotropic 
medium, neglecting the existence of the bottom face of the plate 
except for the cumulative reflection coefficients, to compute the 
delay time, the displacements and the directivity functions 
becomes inappropriate. 

For an SV wave travelling in an isotropic medium, its phase 
velocity is a constant value; that is, its phase velocity is 
directionally independent. However, as shown in Fig. 10, the 
phase velocity of an SV wave travelling in a unidirectional 
fiberglass epoxy composite plate is a function of the angle of 
incidence 6 . Consequently, the directional dependence of the 
phase velocity of an SV wave in the unidirectional fiberglass 
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epoxy composite plate has an effect on the delay time when con- 
ducting acousto-ultrasonic input-output characterization. This 
phenomenon can be observed at low values of dimensionless separa- 
tion 1 / h in Fig. 11, where sharp increases in delay time t n are 
caused by an increase in the number of reflection n from the bot- 
tom face of the plate at a given value of the dimensionless sepa- 
ration £/h. Low values of dimensionless separation f/h correspond 
to an angle of incidence 9 (given in Eq. (49)) whose value is 
less than 60° which is the abscissa of the maximum phase velocity 
C;j_, as shown in Fig. 10. For angles of incidence 9 less than 
60° , the phase velocity increases with increasing angle of 
incidence 9, as shown in Fig. 10. However, an increase in the 
number of reflections n from the bottom face of the plate for a 
given plate thickness h and a given separation distance H results 
in a decrease in the angle of incidence 9 given by Eq. (49), 
thereby increasing the travelling distance according to Eq. 

(50) and decreasing the accompanying phase velocity. The 
increase in the travelling distance Rjj and the decrease in the 
phase velocity account for the sharp increases in the delay 
time t n at low values of dimensionless separation 1/ h. For the 
case of an SV wave travelling in an isotropic medium, the 
increase in the delay time t n caused by an increase in the number 
of reflections n is attributed solely to the increase in the 
travelling distance Rjj. 
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Consider the case of a given number of reflections n. For a 
given plate thickness h, an increase in the separation distance Z 
results in an increases in the travelling distance R^, as shown 
in Fig. 9, and an increase in the angle of incidence 9 given by 
Eq. (49). An increase in the angle of incidence 9 may result in 
an increase or a decrease in the SV-wave phase velocity C^, 
depending on whether the angle of incidence 9 is less than or 
greater than 60°, as shown in Fig. 10. However, as the separa- 
tion distance i approaches infinity, the angle of incidence 9 
will approach 90®, and the phase velocity will reach a limit, 
as shown in Fig. 10. Consequently, (for i/h -* °°) the increase in 
the delay time t n is due almost exclusively to the increase in 
the travelling distance R^, which is similar to the isotropic 
medium case. 

This theoretical investigation provides a step forward in the 
quantitative understanding of acousto-ultrasonic nondestructive 
evaluation (NDE) parameters such as the stress wave factor (SWF) 
in transversely isotropic media. It also provides the potential 
for assisting in the development of more efficient and more 
revealing NDE schemes utilizing wave propagation. 
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z 



Fig. 2 Coordinate system (x,y,z) in the analysis of single 
reflection problem at stress-free plane boundary of 
semi-infinte transversely isotropic medium for angle 
of incidence less than critical angle; z = 0 is plane 
boundary where single reflection occurs, z < 0 is free 
space, and x = 0 is plane of incidence. 
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* Slowness surface for P wave 
; Slowness surface for SV wave 
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Two sheets of slowness surface for P and SV waves in unidirectional fiberglass 
epoxy composite for positive y-z quadrant. 
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Fig. 4 Angle of reflection of reflected P wave versus angle of incidence of 
incident SV wave in unidirectional fiberglass epoxy composite. 


A (sv > 

: Amplitude ratio of reflected SVwave, jj, 



4 ‘ono i 3pnii|dujv 

(d) V (AS) V 


67 


Fig. 5 Amplitude ratios of reflected SV and P waves to incident SV wave versus 

angle of incidence in unidirectional fiberglass epoxy composite for angle of 
incidence less than critical angle. 


P wave 

amplitude decreases 
exponentially with z 



Fig. 6 Existence of reflected surface P wave whose amplitude 

decreases exponentially with z for angle of incidence of 
incident SV wave in unidirectional fiberglass epoxy 
composite, equal to or greater than critical angle 0 C ; 
z = 0 is plane boundary where single reflection occurs, 
z < 0 is free space, and x = 0 is plane of incidence. 
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• Z-component energy flux of ref lected SV wave 
s Z-component energy flux of reflected Pwave 



Angle of incidence, Q\ (degree) 

Energy fluxes normal to plane boundary, assuming radian frequency equal to unity 
for incident SV wave in unidirectional fiberglass epoxy composite for angle of 
incidence less than critical angle. 





» 




_ O O' 00 N- id in 
cvj cvi — “ " 

(S/OJ)j) <0 ‘ 1 D 0 1 3A 9SDL|d 


72 


Fig. 10 Phase velocity versus angle of incidence for SV wave travelling in 
unidirectional fiberglass epoxy composite plate specimen. 
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Fig. 13 Polar diagram for directivity function Dy Z of shear stress Tyz at transmitting 

origin associated with SV wave in unidirectional fiberglass epoxy composite plate 
specimen at frequency of 0.75 MHz for positive y-z quadrant due to applied point 



Polar diagram for directivity function Dy z of shear stress Ty Z at transmitting 
origin associated with SV wave in unidirectional fiberglass epoxy composite plate 
specimen at frequency of 1.50 MHz for positive y-z quadrant due to applied point 
force acting along y-direction. 
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Fig. 16 Polar diagram for directivity function D zz of normal stress x zz at transmitting 
origin associated with SV wave in unidirectional fiberglass epoxy composite 
plate specimen at frequency of 0.75 MHz for positive y-z quadrant due to applied 
point force acting along y-direction. 
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Directivity function D zz along y axis 

Fie. 17 Polar diagram for directivity function of normal stress T zz at transmitting 

or igin associated with SV wave in unidirectional fiberglass epoxy composite plate 
specimen at frequency of 1.50 MHz for positive y-z quadrant due to applied point 
force acting along y-direction. 
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Fig. 19 Polar diagram for directivity function D of shear stress Ty Z at transmitting 

origin associated with SV wave in unidirectional fiberglass epoxy composite plate 
specimen at frequency of 0.75 MHz for positive y-z quadrant due to applied point 
force acting along z-direction. 
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Directivity function D zz along y axis 

Fig. 22 Polar diagram for directivity function D zz of normal stress x zz at transmitting 

origin associated with SV wave in unidirectional fiberglass epoxy composite plate 
specimen at frequency of 0.75 MHz for positive y-z quadrant due to applied point 
force acting along z-direction. 
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